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A new frequency standard has been demonstrated with the aid of a double phase - 
locked loop (TIL) receiver. A superconducting cavity-maser oscillator (SCMO) 
and a hydrogen maser are combined to show the medium-term performance of the 
hydrogen maser together with improved short-term performance made possible by 
the SCMO. The receiver, which generates a 100-MHz signal with reduced noise, 
is phase-locked to ( and may be used in place of) the 100-MIIz hydrogen-maser 
output. The maser signal, 2.69xxx-GHz SCMO output, and a 100-MIIz quartz- 
crystal oscillator are optimally combined by the receiver. A measured two-source 
fractional frequency stability of 2 x 10~ 14 was obtained for a measuring time of 
r = 1 sec, and 1 x 10~ 15 at t = 1,000 sec. The 1-sec value is approximately 10 
times lower than that for hydrogen masers, while the 1,000-sec value is identical 
to hydrogen-maser performance. The design is based on phase-noise models for 
the hydrogen maser, the SCMO, and quartz-crystal oscillators for offset frequencies 
down to 1 x 10- 6 Hz. 


I. SCMO Operation 

The superconducting cavity-maser oscillator (SCMO) 
is an all-cryogenic, helium-cooled oscillator with superior 
stability at short measuring times [1-4]. It differs from 
other superconducting cavity-stabilized oscillator (SCSO) 
designs [5-7] in its use of a very rigid (Q « 10 9 ) sapphire- 
filled stabilizing cavity and in its all-cryogenic design, with 
excitation provided by an ultralow-noise cryogenic ruby 
maser. 

The three-cavity oscillator, which consists of a ruby 
maser, coupling cavity, and a high- Q lead-on-sapphire cav- 
ity, has been discussed previously [2]. Oscillation at a 
frequency of 2.69xxx GIIz results from maser amplifica- 
tion in the ruby material, where a population inversion is 


generated by a 13.1-GIIz pump signal. A match between 
ruby and resonator frequencies is effected by application to 
the ruby of an approximately 500-gauss bias field. Output 
power of the oscillator is > 10 ^ W , more than 1,000 times 
larger than that of the hydrogen maser, which makes pos- 
sible \/l000 « 30 times higher stability at short measuring 
times. 

The SCMO was recently modified to allow its frequency 
to be actively adjusted to match that of the hydrogen 
maser. A coil has been installed on the ruby housing to 
modify the magnetic bias field by application of a DC cur- 
rent. The coil produces a tuning sensitivity of 7 x 10 
per mA, with a range of approximately 10~ 10 without sig- 
nificant heating. This range, although only 1/1000 that 
of a quartz-crystal voltage-controlled oscillator (VCO), 
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is sufficient to accommodate the typical SCMO drift of 
4 x 10~ 13 /day in long-term operation. 


IL Loop Design and Oscillator Noise Models 

A double phase-locked loop (PLL) was designed to op- 
timally combine SCMO and hydrogen-maser stabilities. 
Figure 1 shows a block diagram of that design. Here, 
one loop locks the phase of a quartz-crystal VCO to the 
SCMO, while the second loop locks the SCMO/VCO com- 
bination to a hydrogen maser through the tuning coil in the 
SCMO, Design goals are to preserve SCMO short-term sta- 
bility through the first PLL and optimize the second PLL 
so that long-term performance of the hydrogen maser is 
preserved without significantly degrading SCMO perfor- 
mance at t = 1 sec measuring time. The VCO-SCMO 
loop characteristics, with a bandwidth of approximately 
1 , 000 Hz , do not significantly impact stability performance 
when r > 1 sec. However, the SCM O-hydrogen-maser 
loop does, and so deserves more attention. 

While most measurements of the stability of frequency 
standards are expressed as an Allan deviation <r y {r) in 
the time domain, the results do not directly apply to loop 
design. This is because the stability (for any measuring 
time r) of a source that combines several standards will 
depend on performance of the contributing standards at 
every r. 

However, the performance of any standard can also be 
expressed in the frequency domain as a spectral density 
of phase fluctuations S*(/), and the spectral densities for 
the various standards at any frequency / simply add, with 
multiplicative constants determined by loop characteris- 
tics at that frequency. Roughly speaking, the unity-gain 
frequency for the loop will match the crossover frequency 
between the spectral densities for any two standards being 
combined. 

Although stability measurement data are abundantly 
available for the DSN hydrogen masers, there existed no 
measurement of close-in phase noise (/ < 1 IIz) when the 
authors undertook their receiver design. Thus, they cre- 
ated a phase-noise model for which the calculated Allan 
deviation <r y (r) matches in detail the results of stability 
measurements in JPL’s frequency standards test facility 
(FSTF), 

Details of this model appear in Table 1. Included 
for three separate noise components are! noise type, the 
Fourier frequency window for which the component is dom- 
inant, Allan deviation, and phase spectral density. The 


authors* calculation methodology for generating the Al- 
lan deviation corresponding to any phase-noise model is 
described in the Appendix. 

Figure 2 shows a comparison of this model with the re- 
sults of a very recent calculation of hydrogen-maser close- 
in phase noise using the same raw data from which Allan 
deviations were calculated. These data consist of measure- 
ments of the time for zero-crossings of a 1-Hz difference 
frequency using a 100-MHz -f 1-Hz reference. The differ- 
ence between the model and the measured data is within 
3 dB for the region from 1 x 10~ 4 to 1 x 10 -1 . 

Measurements of the spectral density of phase fluctu- 
ations were available for the other two frequency sources. 
The authors had previously made measurements of the 
phase noise for the SCMO/SCSO combination for offset 
frequencies down to 0.01 Hz, and in the absence of any bet- 
ter information, they ascribed half of the measured noise 
to the SCMO. Other tests of the SCMO with a hydro- 
gen maser as reference indicate a somewhat lower value 
for offsets below 0.01 Hz. Information for the lOO-MIIz 
quartz-crystal oscillator was obtained from the manufac- 
turer. 

Figure 3 shows a plot of the phase noises, measured 
and inferred, for the three oscillators. A smaller contri- 
bution due to input noise for the operational amplifier in 
loop no. 2 is also shown. Crossover frequencies of approx- 
imately 0.04 Hz and approximately 1,000 Hz are readily 
identifiable. 


HL Loop-Parameter Optimization 

As shown by the dotted line in Fig. 3, calculated per- 
formance for the combined source closely approximates 
the best of the three sources at every frequency. How- 
ever, deviations from “ideal” performance are significant, 
as Fig. 4 shows in expanded form, for second-order loops 
with various loop bandwjdths and damping factors. Pa- 
rameter optimization requires a measure of “goodness” for 
loop action. It is often the case that performance is op- 
timized for a time period which is much longer than that 
associated with action of the loop. In this case such an 
optimization measure can be, for example, rms deviation 
of the phase error of the following oscillator from the one 
being followed. In the case at hand, such a single high- or 
low-frequency measure is not useful because frequencies of 
concern lie both above and below the crossover frequency. 

As shown in Fig. 4, performance deviates significantly 
from the ideal for frequencies above and below the cross- 
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over frequency itself. Some applications of the com- 
bined source will, for example, emphasize high- over low- 
frequency performance (or short-term stability for a com- 
bined source at the expense of long-term stability). The 
authors have chosen instead a balanced optimization crite- 
rion, which both equalizes and minimizes peak deviations 
on either side of the crossover. 


Successive steps approaching this optimized condition 
are shown in Fig. 4. A first guess— matching the loop 
bandwidth to the crossover frequency (0.04 Hz) with a 
damping factor of l/\/2 — resulted in a large low-frequency 
peak, substantially worse than any shown in the figure. 
Equivalent high- and low-frequency performances can be 
observed in the examples on a diagonal from upper right to 
lower left. Optimized performance, as shown at the lower 
left, results primarily from an increased damping factor. 

The authors’ calculations show continuing small im- 
provement (< 0.2 dB) as the damping factor is further 
increased. This indicates that on the basis of this cri- 
terion alone, an infinite damping factor, corresponding 
to a first-order loop, would be preferable. However, the 
larger frequency drift rate of the SCMO, as compared with 
the hydrogen- maser, necessitates the use of a second-order 
loop to prevent a slow buildup of phase error, which cor- 
responds to a frequency offset, between the sources. The 
value of 3 for the damping factor is chosen to give good 
noise performance and drift immunity combined with rel- 
atively rapid loop response. 

While spectral densities give a complete representation 
of performance of the standard, performance in the time 
domain is a more familiar and widely used measure of per- 
formance. Figure 5 shows a calculated Allan deviation for 
the composite spectral density for the combined sources. 
It is surprising that stability near the crossover point here 
is a little better than that of either source, while it closely 
follows the SCMO and the hydrogen maser at short and 
longer times, respectively. The slight improvement at the 
crossover occurs because the Allan deviation for either 
component standard depends significantly on frequencies 
both above and below the crossover frequency, and the 
composite has the advantages of both at the crossover it- 
self. 

The calculated combined source shows a stability of 
1 x 10 -14 at 1 sec and 1 x 10 -15 from 1,000 sec onward. 
The Allan deviation of the combined source follows very 
closely the best source in the time window of 1-10,000 sec. 


IV. Measurements 

To confirm the stability of the combined source, a sec- 
ond identical system will be needed. Presently available 
are hydrogen-masers to test performance at longer measur- 
ing times and the SCSO for short-term characterization 
[5]. Figure 6 shows measurements that demonstrate the 
short-term stability of the combined source at a test fre- 
quency of 100 MHz using the SCSO as a reference. During 
this test, the poor long-term stability of the SCSO refer- 
ence substantially degraded the measured values for times 
r > 30 sec. However, results at shorter times are as much 
as 10 times lower than for two hydrogen masers in a similar 
test. The slight increase for r < 3 sec has previously been 
identified as characteristic of the SCSO reference [6], which 
indicates performance for the authors’ combined source of 
<Ty(r) < 1 x 10 -14 at t - 1 sec. 

Figure 7 shows the Allan deviation of fractional fre- 
quency fluctuations for the SCMO locked to a hydrogen 
maser, using a second hydrogen maser as reference. A 
special feature is the stability of 6 X 10 14 at 1 sec, a value 
that is one maser’s stability. The transition from “one- 
maser” noise below t = 30 sec and “two-maser noise for 
longer times causes the slight kink shown in the figure. 

V. Conclusions 

The authors have designed, built, and demonstrated a 
receiver that enables the SCMO to operate in conjunc- 
tion with a hydrogen maser to form a new standard, one 
which combines their complementary performance over the 
range of measuring times from 1 to 10,000 sec. Calcu- 
lated performance for the combined source shows spectral 
performance 5^(/) within 1.4 dB of the best of the two 
sources over the frequency range 10 6 Hz < / < 10 Hz. 
Allen deviation is calculated to be within 7 percent of the 
best of the two sources for measuring times of 1 sec < 
r < 10 4 sec, and somewhat surprisingly, slightly better 
than either standard at their crossover. Performance for 
the combined standard was demonstrated in separate ex- 
periments using different ultrastable frequency sources to 
be at least as good as 2 x 10“ 14 at 1-sec measuring times 
and 1 x 10" 15 at 1,000 sec. Significant aspects of this ex- 
periment are a new time window for scientific experiments 
and a unique demonstration that combines two different 
types of ultrastable microwave oscillators. 

The combined SCMO/hydrogen-maser source under- 
went a field test at Goldstone Deep Space Communica- 
tions Complex for 72 days (May-July 1991). The perfor- 
mance was the same as reported above. This exercise was 
to prepare the SCMO for a gravitational wave detection 
experiment in connection with Galileo in 1992. 
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Table 1 . Contributions to the total noise of the DSN hydrogen maser used for loop design 


Noise type 

Fourier frequency window 

Allan deviation 

Phase spectral density 

Flicker frequency 
White frequency 
Flicker phase 

1 x 10“ 6 to 4 X 10~ 4 
4 X 10“ 4 to 4.5 X 10“ 2 
4.5 X 10" 2 to 1 X 10 _1 

1 X 10- 15 

3 X 10 - 1 Vt 1/2 

2 X 10 _13 /t 

7.5 X 10“ 15 // 3 
1.8 X 10 - 11 // 2 
4 X 10 — 10 // 



Fig. 1. Schematic of SCMO double-loop receiver at 100 MHz. The two loops are VCO/SCMO and SCMO/hydrogen 
maser. The required input is the 100-MHz hydrogen-maser signal, and the output has four different frequencies- 
100, 10, 5, and 0.1 MHz. M 
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Fig. 2. A close-ln phase-noise plot. The solid line shows the model 
maser; two sets of measured pair-data were shown. The difference 

region from 1 xIO -4 to 1 xIO 1 Hz. 


of a single hydrogen 
Is within 3 dB for the 



Fig. 3. Phase-noise plot of the three oscillators at 100 MHz. Calculated receiver 
shown by the dotted line, closely follows the lowest noise of the three oscillators In the Four 
frequency window of IxIlT 5 to 10 4 Hz. Noise due to the operational amplifier in loop no. 2 
is shown to be much lower than oscillator noise contributions. 
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Fig. 4. Optimization of the loop parameters Is accomplished by minimizing calculated excess 
phase noise at the receiver output compared with the best of the constituent oscillators at 
any given offset frequency, Bandwidths for the various graphs Increase from left to right 
and damping factor increases from top to bottom. Horizontal and vertical scales are Identical 
for each subgraph. Best overall performance, as shown at the lower left, corresponds to a 
bandwidth of 0.04 Hz (a value equal to the SCMO-hydrogen-maser crossover frequency) and 
a relatively large damping factor. 
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Fig. 5. Allan deviation derived from combined phase noise data 
(Fig. 3). The expected performance of the combined signal is 
1 x 10” 14 at 1 sec and 1 x 10“ 15 at 1,000 sec and beyond. 
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Fig. 7. Two sample Allan deviations of the SCMO locked to one 
hydrogen maser; the second maser is used as a reference. The 
vertical line shows the 30-sec loop time constant. The transition 
from "one-maser" noise below r — 30 sec and "two-maser" 
noise for longer times causes the slight kink shown in the figure. 


Fig. 6. Two sample Allan deviations of the SCMO tested with 
SCSO at 100 MHz. The measured frequency stability is 2 x 10“ 
at 1 sec. Structure below r = 100 sec Is primarily due to the 
SCSO reference. 
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Appendix 

Calculation of the Allan Deviation 


The Allan deviation of variation of fractional frequency 
cr^(r) for any radio frequency (RF) signal can be derived 
from its spectral density of phase noise S<p(f) by using [8] 


2 

y 



oo 


C(r) x S#(f) x sin 4 (?r rf)df 


typically — 61n(10) and 4 ln( 1 0) , which corresponds to an 
integration over frequencies from 10 -6 Hz to 10 4 Hz, 


Elimination of rapid variation is accomplished by re- 
placing 


sin 4 (7rr/) 


where 


C(r) = 


2 

7 r 2 r 2 z/2 


by 


Die 



x 




and Vo is the RF frequency. In order to use the numerical 
integration routines in a popular mathematics program [9] 
to evaluate this relatively straightforward integral, it is 
necessary to “condition” the problem in two ways. 

First, because significant contributions to the integral 
are due to frequency components that vary by many orders 
of magnitude, a change in variable is required, which re- 
sults in an effectively logarithmic frequency scale. Second, 
calculation at high frequencies, as compared with 1/r, is 
hampered by rapid variation in the sin 4 term, which should 
be replaced by its average value (3/8) as the frequency is 
raised. It is important for the function that does the re- 
placement to show very complete elimination of the aver- 
age term as the argument nrf goes to zero so as not to 
dominate the (7rr/) 4 dependence of the sin 4 term. 

The change of variable required is straightforward, re- 
placing the frequency / with e 3 and the differential df 
with e 9 dg. Lower and upper limits to the integration were 


where n represents the number of cycles of variation in the 
sin 4 term that are to be explicitly integrated, and Dk is a 
decay function defined by 

Dk(x) = e~ x]n W 

for x > 1, and 

_(\nLl 1) 

Dk(x) = l-e V * ) 

for x < 1. Features of this decay function are extremely 
rapid approach to 1 for small x , rapid approach to 0 for 
large x, and continuous first derivative at the crossover 
point x = 1. The authors found n = 3 cycles sufficient to 
reduce fractional errors to less than 10“ 4 . 

Taken together, these substitutions allow a difficult in- 
tegral to be accurately evaluated by means of “canned” 
integration routines. 
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